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ABSTRACT 

Thermal instability in an electron-ion magnetized plasma which is relevant in the 
intragalactic medium (IGM) of galaxy clusters, solar corona, and other two-component 
plasma objects is investigated. We apply the multicomponent plasma approach when the 
dynamics of all species is considered separately through the electric field perturbations. 
General expressions for the dynamical variables obtained in this paper can be applied for 
Jg j a wid e rang e o f a8tr0 p hy8ical a„ d laborat0ry Plasmas a ls0 c„n, ining — and du8t 
grains. We assume that background temperatures of electrons and ions are different and 
include the energy exchange in the thermal equations for the electrons and ions along with 
the collisional momentum exchange in the equations of motion. We take into account 
the dependence of collision frequency on the density and temperature perturbations. The 
cooling-heating functions are taken for both electrons and ions. A condensation mode 
of thermal instability has been studied in the fast sound speed limit. A new dispersion 
relation including the different electron and ion cooling-heating functions and other effects 
mentioned above has been derived and its simple solutions for growth rates in the limiting 
cases have been found. We have shown that the perturbations have an electromagnetic 
nature. The crucial role of the electric field perturbation along the background magnetic 
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field in the fast sound speed limit has been demonstrated. We have found that at conditions 
under consideration, the condensation must occur along the magnetic field while the 
transverse scale sizes can be both larger and smaller than the longitudinal ones. The 
results obtained can be useful for interpretation of observations of dense cold regions in 
astrophysical objects. 

Key words: conduction - galaxies: clusters: general - instabilities - magnetic fields - 
plasmas -waves 
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1. INTRODUCTION 

If a medium in a thermal equilibrium can become cooler due to radiation and fluid 
contraction, it can get unstable, leading to formation of density condensations with 
lower temperature than in the surrounding medium (Parker 1953; Field 1965). This 
instability called the thermal or radiation-condensation one has been studied for more 
than five decades for astrophysical objects and plasma physics applications (for reviews 
see, e.g., Meerson 1996; Vazquez- Semadeni et al. 2003; Elmegreen k Scalo 2004; Cox 
2005; Heiles k Crutcher 2005). Many papers in astrophysical literature considered the 
thermal instability in the neutral (Field 1965; Begelman k McKee 1990; Hennebelle k 
Perault 1999; Koyama k Inutsuka 2000, 2002; Burkert k Lin 2000; Kritsuk k Norman 
2002; Sanchez-Salcedo et al. 2002; Audit k Hennebelle 2005; Vazquez- Semadeni et al. 
2006; Hennebelle k Audit 2007 and references given above) and magnetized interstellar 
medium (ISM; Field 1965; Hennebelle k Perault 2000; Stiele et al. 2006; Fukue k Kamaya 
2007; Inoue k Inutsuka 2008; Shadmehri et al. 2010), solar corona, where prominences 
are formed (e.g., Field 1965; Nakagawa, 1970; Heyvaerts 1974; Mason k Bessey 1983; 
Karpen et al. 1988), planetary nebulae (e.g., Field 1965), galaxy clusters and intragalactic 
medium (IGM; Field 1965; Mathews k Bregman 1978; Balbus k Soker 1989; Loewenstein 
1990; Balbus 1991; Bogdanovic et al. 2009; Parrish et al. 2009; Sharma et al. 2010). 
This instability was studied for rotating (e.g., Field 1965; Nipoti 2010), expanding (e.g., 
Field 1965; Gomez-Pelaez k Moreno-Insertis 2002), and dynamical systems with steady 
or time- dependent flows and nonstationary background parameters (e.g., Mathews k 
Bregman 1978; Balbus 1986; Balbus k Soker 1989; Burkert k Lin 2000). The nonlinear 
stage of the thermal instability resulting in formation of nonlinear structures (localized 
clouds) in the ISM and solar prominences was investigated in (e.g., Trevisan k Ibahez 2000; 
Sanchez-Salcedo et al. 2002; Yatou k Toh 2009) and (Mason k Bessey 1983; Karpen et al. 
1989; Trevisan k Ibanez 2000), respectively. 
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In majority of papers studying the thermal instability of astrophysical objects in the 
magnetic field, one uses the one-fluid ideal magnetohydrodynamics (MHD). In some papers, 
the two- fluid model with the ideal magnetic induction equation has been treated (e.g., 
Fukue & Kamaya 2007; Inoue & Inutsuka 2008). The non-ideal effects in the induction 
equation have been included in (e.g., Heyvaerts 1974; Stiele et al. 2006; Shadmehri et al. 
2010). 

In the plasma physics literature devoted to the thermal instability, one considers besides 
the one-fluid (e.g., Bora & Taiwar 1993) and two-fluid (e.g., Birk 2000) ideal MHD (from 
the point of view of the form of the magnetic induction equation) also a multicomponent 
media where the presence of electrons, ions, dust grains, and neutrals is taken into account 
(Kopp et al. 1997; Birk & Wiechen 2001; Pandey & Krishan 2001; Pandey et al. 2003; 
Shukla & Sandberg 2003; Kopp & Shchekinov 2007). Such effect as dust charge variation 
is also included when studying the thermal instability (ibanez & Shchekinov 2002; Pandey 
& Krishan 2001; Pandey et al. 2003). Analytical investigation of thermal instability in the 
multicomponent magnetized media with such physical effects as collisions between different 
species, ionization and recombination, dust charge dynamics, gravity, self-gravity, rotation, 
and so on is a very challenging problem. General basic equations have a very complex form 
(see, e.g., Kopp et al. 1997). Therefore, one usually treats simplified models considering, for 
example, potential perturbations in nonmagnetized (Kopp et al. 1997; Pandey & Krishan 
2001; Ibanez & Shchekinov 2002; Pandey et al. 2003; Shukla & Sandberg 2003; Kopp k 
Shchekinov 2007) and magnetized (Kopp et al. 1997; Shukla & Sandberg 2003) plasmas. 

Neglecting the energy exchange between species in the thermal equations is one of the 
usual simplifying assumptions. It is justified when a collisional coupling is weak or when 
the energy exchange frequency is sufficiently large, so that the temperatures of species are 
equal. In a general case, one can obtain considerable complications. Nevertheless, some 
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authors took into account the energy exchange for thermal instability in the two-fluid MHD 
framework (e.g., Birk 2000; Birk & Wiechen 2001). 

The terms describing the energy exchange contribute to the dispersion relation not 
only through the temperature (or pressure) perturbations (Birk 2000; Birk & Wiechen 
2001) but also through the perturbation of collision frequency which depends on the density 
and temperature. This effect is important when background temperatures of species are 
different. Such a situation can occur, for example, in galaxy clusters (see Markevitch et al. 
1996; Fox & Loeb 1997; Ettori & Fabian 1998; Takizawa 1998). 

In the present paper, we investigate the thermal instability in the electron-ion 
magnetized plasma which is relevant in the IGM of galaxy clusters, solar corona, and other 
two-component plasma objects. We apply the multicomponent plasma approach when the 
dynamics of all species is considered separately through the electric field perturbations (the 
E-approach; see, e.g., Nekrasov 2009 a, 2009 b, 2009 c; Nekrasov & Shadmehri 2010, 2011). 
General expressions obtained in this paper can be applied for a wide range of astrophysical 
and laboratory plasmas. We assume that the background electron and ion temperatures are 
different and include the energy exchange in the thermal equations for electrons and ions. 
We take into account the dependence of collision frequency on the density and temperature 
perturbations. Different cooling-heating functions are assumed for electrons and ions. We 
include neither ionization and recombination effects nor gravity. Some expressions for 
electron and ion perturbations are obtained in the general form which can be used for 
other species (dust grains and neutrals). Here, we treat the condensation mode of thermal 
instability in the fast sound speed limit. We derive the general dispersion relation taking 
into account the effects mentioned above and find its simple solutions for the growth rates 
in the limiting cases. 

The paper is organized in the following manner. In section 2, we give the fundamental 
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equations used in this paper. The equilibrium state is considered in section 3. Equations for 
temperature perturbations are obtained in section 4. In section 5, we give specific conditions 
for further consideration. In section 6, equations for components of velocity perturbations 
are given in the fast sound limit. Components of the perturbed current are calculated in 
section 7. These components for the simplified collision contribution are given in section 8. 
In section 9, we derive the dispersion relation. Its limiting cases are considered in section 
10. We discuss the obtained results in section 11. The possible astrophysical implications 
are considered in section 12. Summing up of the main points is given in section 13. 



2. BASIC EQUATIONS 



The fundamental equations that we consider here are the following: 



8V, +vrW, = -^- + P, + JLv < xB, (1) 



dt m jUj rrijC 



the equation of motion, 

^ + V • n jVj = 0, (2) 

the continuity equation, 

dT- 1 

-± + v, • VT, + ( 7 - 1) T,V ■ v, = - ( 7 - 1) —Li (m, Ti) + v\ e (n e , T e ) (T e - T t ) (3) 
dt rti 

and 

c)T 1 X 

-T^+Ve-VT e +(7 - 1) T e V-v e = - (7 - 1) —V • q e -( 7 - 1) — C e (n e , T e )-v e ei (m, T e ) (T e - T,) 

(4) 

are the temperature equations for ions and electrons. In Equations (1) and (2), the index 
j — i,e denotes the ions and electrons, respectively. The force Fj in Equation (1) is given 
by 
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E - (Vi - V e ) , 



(5) 




mi 
m e 



E - Z/ ei (v e - Vi) . 



Other notations in Equations (l)-(5) are the following: qj and rrij are the charge and mass 
of species j — i,e, Vj is the hydrodynamic velocity, rij is the number density, pj = rijTj 
is the thermal pressure, Tj is the temperature, u ie (u ei ) is the collision frequency of ions 
(electrons) with electrons (ions), vf e {n e ,T e ) = 2u ie {y e ei (rii, T e )) is the frequency of the 
thermal energy exchange between ions (electrons) and electrons (ions) (Braginskii 1965), 
niU? e (n e ,T e ) = n e v e ei (n i: T e ), 7 is the ratio of specific heats, E and B are the electric and 
magnetic fields, and c is the speed of light in vacuum. The value q e in Equation (4) is the 
electron heat flux (Braginskii 1965). As for the latter, we will consider a weakly collisional 
plasma when the electron Larmor radius is much smaller than the electron collisional mean 
free path. In this case, the electron thermal flux is mainly directed along the magnetic field, 



where Xe is the electron thermal conductivity coefficient and b = B/5 is the unit vector 
along the magnetic field. In other respects, a relation between cyclotron and collision 
frequencies of species stays arbitrary in general expressions considered below. We only 
take into account the electron thermal flux (6) because the corresponding ion thermal 
conductivity is considerably smaller (Braginskii 1965). We also assume that the thermal 
flux in the equilibrium is absent. The cooling and heating of plasma species in Equations 
(3) and (4) are described by function £j(rij,Tj) = n^Aj (Tj) — rijTj, where A,- and Tj 
are the cooling and heating functions, respectively. The form of this function has some 
difference from the usually used cooling- heating function £ (Field 1965). Both functions 
are connected with each other via equality Cj (rij,Tj) = mjTij£j. Our choice is analogous 



q e = - 



Xeb(b-V)T £ 



ei 



(6) 
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to that as in (Begelman & Zweibel 1994; Pandey & Krishan 2001; Shukla & Sandberg 
2003; Bogdanovic et al. 2009; Parrish et al. 2009). The function A,- (7}) can be found, for 
example, in (Tozzi & Norman 2001). 

Electromagnetic equations are Faraday's 



and Ampere's 



V x B =-j (8) 

c 



laws, where j = J2j 9j n i v i- We consider wave processes with typical timescales much larger 
than the time the light spends to cover the wavelength of perturbations. In this case, 
one can neglect the displacement current in Equation (8) that results in quasi-neutrality 
for both the electromagnetic and purely electrostatic perturbations. The magnetic field 
B includes the background magnetic field B , the magnetic field B 0c «r of the background 
electric current (when it presents), and the perturbed magnetic field. 

For generality, we assume in the meanwhile that n« and n e are different, having in mind 
that some solutions obtained below can be applied for multicomponent plasmas. 



3. EQUILIBRIUM STATE 

At first, we will consider an equilibrium state. We assume that the background 
velocities of species are absent. Here, we do not involve an equilibrium inhomogeneity. 
Then, the thermal equations (3) and (4) in the equilibrium take the form 

(7 - 1) — A (n i0 ,T i0 ) - z4(n e0 , T e0 ) (T e0 - T i0 ) = 0, (9) 

(7 - 1) — C e (n e0 , T e0 ) + v £ ei (n i0 , T e0 ) (T e0 - T i0 ) = 0. 



4. LINEAR EQUATIONS FOR THE TEMPERATURE PERTURBATIONS 

We now consider Equations (3) and (4) in the linear approximations. Using Equations 
(2) and (9) for ions, we find 



D u T a - D 2i T el = C U V ■ v a - C 2l V ■ v e i, 



(10) 



where the index 1 denotes the perturbed values. The following notations are here 
introduced: 



9 n n V 

dt + ilTl + ate ) w 

d 



D 2 i — {&Tie + &ie) q^, 



C\i — T i0 

C 2 i = (T e Q — Tjo 



-(7-1)1 + ^-2^^1* 

Cft liO 



where the index denotes the unperturbed values. Analogously, we obtain for the electrons 



D le T el - D 2e Tn = Ci e V • v ei + C 2e V • v a , 



(12) 



where 



d_ 

dt 



Die — \ — + + ^Te + ^Tei + ^ei 



d_ 



(13) 



D 2e — VL e i 



d_ 



Cle — T e o 
C 2e = fl e i (T e Q — Tjo) . 



(7 - !) ^7 + ^ne + = ^ei 

Ot J eO 
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In notations (11) and (13), we have introduced the following frequencies: 

^-fr- 1 '^' (14) 

dC e (n e0 ,T e0 ) dCi(n i0 ,T i0 ) 
QTe = (7 " 1} n e0 dT e0 ' = (7 " 1} n^T 40 ' 

n~ = (7 - 1) ^; (n : 0;Teo) , = (7 - 1) ^; (n :°' T * o) , 

&ei = v £ ei (Uio, T e0 ) , fii e = z/f e (n e0 , T e0 ) , 

_ <9zA(n t0 ,T e0 ) disf e (n e0 ,T e0 ) 

^ L Tei — ^ UeO ~ 1 iO) , "Tie — ^ U eO — 1 iO) ■ 

CI e Ol e o 

We assume that the background magnetic field B is directed along the z-axis. In notations 
(11) and (13), we have used an equilibrium state and the fact that v e ei (n i0 ,T e0 ) ~ n^o 
and vf e (n e o,T e0 ) ~ n e0 . We see from Equations (10) and (12) that the temperature 
perturbations are connected with a velocity divergence. Solutions for T el and T a are given 
by 

DT el = G 1 V • v e i + G 2 V • va, (15) 
.DTii = G 3 V ■ v e i + G 4 V • Vji, (16) 
where the following notations are introduced: 

D = D li D u -D 2i D 2e , (17) 

G\ = DnCi e — D2 e C2i, 
G2 = D\iC2e + D2 e C\i, 

G3 — D 2 iCi e — D\ e C2ii 

G4 = D\ e C\i + D2iC2e- 



To find the temperature perturbation Tji, we have to calculate expressions for V • Vji. 
General equations for the velocity v^i and V • Vji are derived in the Appendix, where 
expressions for D and Gi, 1 — 1,2, 3, 4, are also given. In their general form, the components 
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of Vji are very complex. Therefore to proceed further analytically, here we restrict ourselves 
to a limiting case in which the dynamical frequency d/dt is smaller than the sound 
frequency Thus, we will consider sufficiently short-wavelength perturbations along the 
magnetic field (see below). Some additional simplifying conditions which are satisfied in 
magnetized plasmas are also used. 



5. SPECIFIC CASE: FAST SOUND SPEED LIMIT 

Equations (A27), (A28) and (A30) are written in their general form which allows us to 
consider different simplified specific cases corresponding to the real astrophysical conditions. 
We further consider the case in which 



dz J \ dy 2 dz 2 J dt 2 

Since usually io 2 ci ^> d 2 /dt 2 in magnetized plasmas, condition (18) is satisfied for a wide 
range of the transverse wavelengths of perturbations (d/dy)" 1 in comparison with the 
longitudinal wavelengths (d/dz)" 1 . The other simplifying condition is 

, [T e0 T t 4 t + (T e0 + T l0 ) 2 n ie ] fd 2 d 2 \ 



m t [T e0 f t + (T e0 + T t0 )Q te ] \dy 2 dz 2 
This inequality means that the effective ion Larmor radius is much smaller than the 
perturbations wavelengths. 

Below, we consider the case in which the perturbation frequency is much smaller than 
the sound frequency, i.e., 

f_ [T e0 T t0 | + (T e0 + T t0 ) 2 a e ] d 2 
dt 2 mi [ Te0 f t + (T e0 + T i0 )n ie ] dz 2 ' 

This condition can be written in the form d 2 /dt 2 c 2 d 2 /dz 2 , where c s can be considered 
as the effective sound velocity. Inequalities (19) and (20) are given in the approximate form 
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to unite two cases, > (<)Q ie)ei . We note that the case (18) is followed from conditions 
(19) and (20). 

Under conditions (19) and (20), we find, using notations (A10), (A30), and (A31), the 
following equations for P i el (see Equations (A27) and (A28)): 



d 2 Pg = / 1 &Q ilv | d 2 F ilz \ | / 1 d s Q ely | d 2 F elz 

dz 2 1 \tUci dydt 2 dzdt J 6 \co 2 e dydt 2 dzdt 

d 2 Pei / 1 d*Q ely d 2 F elz \ ( 1 d s Q ny d 2 F ilz 

dz 2 " 1 e) \co 2 ce dydt 2 + dzdt J +Pl \u 2 ci dydt 2 + dzdt ' ' 



where 



rrii 



d\ ( <9 V 2 d 2 



Qj = _l, OTe0 - Gl -l - ^, (23) 



Tn ( d \ ( d ' \ 2 

m e / 9 \ " 2 d 3 

-Crq 



K" 3 \dzJ dt r 

Below, we will need the correction terms proportional to a e>i <C 1 and (3 e;i <C 1. We 
note that the contribution of the correction term proportional to Q e i y in Equation (21) is 
at least of the order of m e /rrii in comparison with the term proportional to Qa y . However 
for convenience, we keep this small term for the symmetry with equation for P el . The ratio 
of contributions of the correction terms proportional to F e \ z and F iiz in Equation (21) is of 
the order of 

(T e0 + t i0 ) n ie 

TeOgi + (TeO + ^io) 

The same for terms proportional to Qn y and Q e \ y in Equation (22) is 

nii (T e0 + T i0 ) Vt ei 
m e Tio-TTT + (T e0 + T i0 ) VL ei 
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and for the ratio of terms ~ F ilz and F e i z , we have 

(T e0 + r i0 ) n 

e 



6. EQUATIONS FOR COMPONENTS OF VELOCITIES v; el 

We now obtain equations for components of velocities Vj je i, using Equations (21) and 



(22). 



6.1. Equations for v i: 



i,ely 



From Equations (A3), (21), and (22), we find, using notations (A6), 



Vily = 



to, 



-F, + ± 

± ilx \ 3 



d 2 f d 



dy 2 \dz 



ooz. 



d 2 fd 



-2 



+ A 



l 



dy 2 \dz 



d\~ 2 d 5 F f 



d 2 F tla 
dt 2 

1 



1 OF,, 



ily 



u 2 dt 



d 3 F tly 1 fd 

dt 3 e u ce Uci \9z 



"2 a4 



d F e i x 
dy 2 dt 2 



ely 



+ 



d_ 

^ \ dz J d v 2dt3 u li \ dz 

i f d \ _1 <9 2 

(ajFji 2 + /3 e F e i 2 ) , 



- 1 ^2 



flails 

dydt 



u 2 ci \dz 



dydt 



1 1 



d 2 f d x 1 



d 2 F elx 1 9F el2/ 



<% 2 



^ fd 



dy 2 \dz 

2 



<9 3 F 



ely 





' d 


Tl 

ce 


ydz 




-1 




) 



at 3 

d 5 F ily 1 



dz 



- 2 ^4 



•9 Fjix 
dy 2 dt 2 



j 2 

ce 



d 

dz 



-1 x2 



<9 2 F el2 
dydt 



dydt 



{a e F elz + piFiu) 



(24) 



(25) 



We see that these equations are obtained one from another by changing i -H- e. We note 
that aj, j3j ~ mj. However, we keep in these equations some small terms for the symmetry 
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of both equations. The terms proportional to are needed in equations for v^ e i x . We 
also note that a, t ~ {mi/T ) (d/dz)~ 2 d 2 /dt 2 and & ~ a^/ (d/dt + ft ei ) , if T e0 ~ T i0 ~ T . 
For these estimations, we have used expressions (A18), (A25), (A26), and (A29), assuming 



6.2. Equations for i^eis 

Equations for v i:Slx are found from Equation (A2) by using Equations (24) and (25). 



6.3. Equations for v^ elz 

From Equations (A7), (21), and (22), we obtain, using notations (A6), 

d 2 v ilz 1 d 2 F ilx d 2 ( on „ 8 e „ \ . . 

•/;,,. • —/•,;, (26) 



u 2 ci dydt 2 dydt 2 \u 2 ci lly u 2 ce " 



9 

+ — (aiFiia + (3 e F elz ) , 



(9 2 w e i 2 1 d 2 F e i x d 2 ( a e „ ft ^ , 

— Felx + — i^lJ (27) 



1 <9 3 F el , , «9 3 ( a e „ , A 



+ — (a e F elz + piF ilz ) . 
oz 

We note that in Equation (26) ((27)) the terms proportional to (3 e F e i x , y (o£ e F e i x , y ) are small 
as compared with a.iFn x ^ (f3iFn Xty ). However, we also keep them for the symmetry of these 
equations. 
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7. COMPONENTS OF CURRENT 

We now find components of the linear current ji = ^^<?j^jo v ji- It is convenient to 

j 

calculate the value 4ir(d/dt) 1 ji. We further consider the electron-ion plasma in which 
n e o = nio, Qe = —Qi- In our calculations, we will use an equality m e u ei = miV ie . From 
Equations (A2), (24)- (27), and (5) in the linear approximation, we find 



d 

4vr ( — ) ji x = a xx E lx — a xy E ly + a xz E\ z 



b xx (v ilx - v elx ) + b xy (vny - v ely ) - b xz (v ilz - v elz ) , 



(28) 



-l 



jly 0"yxE\ X ~\- QyyE\y dy Z E\ z 

~ by X (V ilx ~ V elx ) - byy (VHy ~ V A y) + by Z (v ilg ~ Uglz) , 



(29) 



d 

4vr ( — j ji z — —a zx E\ x — a zy E\ y + a zz E\ z 

+ b zx (v ilx - v elx ) + b zy (vny - v e i y ) - b zz (vn z - v elz ) . 



(30) 



Here, the following notations are introduced: 



UJ, 



i + 



d 2 fd 



dy 2 \dz 



-2 



UJ, 



pi 

i &xy 0"yx 3" 
UJ™ 



1 + 



d 2 fd 



dy 2 \dz 



-2 



d_ 

di' 



UJ, 



a *-^-)dy-{dzm) ' a 



UJ 



pi 

yy — y , ">yz ~ ^zy 
UJ„; 



* ~ ^ dy 



UJ, 



{dzdt) '° 2 



i Q>yz ^z 
2 



d fd^ 1 

j 2 

CI 



uj 2 ci dy \dz 



d 



( t^- j Mi = ai~u ie , 



1/9 

where uj. p i = (Airriioq 2 / rrii) is the ion plasma frequency. The value 5 in expression for a z 
has the form 

(32) 



(31) 



rrii 
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8. SIMPLIFICATION OF COLLISION CONTRIBUTION 

The relationship between u ce and u ei or cu ci and u ie (that is the same) can be arbitrary 
in Equations (28)- (30) (except of that in the thermal conduction). We further proceed by 
taking into account that d/dt <C u ci . In this case, we can neglect the collisional terms 
proportional to b xy and b yx (see notations (31)). However, a system of Equations (28)-(30) 
remains sufficiently complex to find ji through Ei. Therefore, we further consider the case 
in which the following condition is satisfied: 



1 > 



v iP d ( d 2 d 2 
+ 



d 



(33) 



bj 2 ci dt \ dy 2 dz 2 J \dz 
It is clear that condition (33) can easily be realized. In this case, we may neglect collisional 
terms proportional to b XjZX and b y>zy . Then, the system of Equations (28)-(30) takes the 
form 

d^ 1 



Here, 



47T 



jlx £xxE\ x &xyE\ y ~\- E XZ E\ Z) 



4^ ( J jly Ey X E\ X + EyyE\y Ey Z E\ z ^ 



47T 



jlz £zxE\ x £zyE\ y + E ZZ E\ Z . 



(34) 
(35) 
(36) 



&xx &xxi ^-xy 



^xy 



^ie a X z d 2 ( d_ 
u%(l + d) dydt \dz 



£yx 



Q-yx / 2 



Vie O'zx 



d 2 { d 



ur ci (1 + d) dydt \dz 



' e yy ~ a yy e yz 



(1 + dY 



^yz 



z y 



(i + d)' y (i + dy zz (i + dy 



(37) 



where 



d = a--—*-— = v ip m~b — 



"zz o Oj. 

u 2 pi dt 



(38) 



-17- 



Parameter d defines the collisionless, d <C 1, and collisional, d >> 1, regimes. Below, we will 
derive the dispersion relation. 

9. DISPERSION RELATION AND ELECTRIC FIELD POLARIZATION 

We will further consider Equations (34)- (36) in the Fourier-representation, assuming 
that perturbations have the form exp (ikr— iuit). Then using Equations (7) and (8), we 
obtain the following system of equations: 

(n 2 — e xx ) E lxk + e xy E lyk — e xz E lzk = 0, (39) 
—tyxEixk + [n 2 z — Eyy) E lyk + (—n y n z + e yz ) E lzk = 0, 
£ Z xEi xk + (—n y n z + e zy ) E lyk + (n 2 y — e zz ) E lzk = 0, 

where E^ is the Fourier-image of the electric field perturbation, n = kc/w. The index k by 
Eifc is equal to k = {k,u}. For the Fourier-images of operators and d, we keep the same 
notations. In general, we see that the longitudinal electric field E Xzk ~ E Xx ^ yk inevitably 
arises when k y ^ and n 2 y — e zz ^ 0. The dispersion relation can be found by setting the 
determinant of the system (39) equal to zero. Contribution of terms proportional to e xy and 
e yx into dispersion relation is small as compared with that of terms ~ e xx ,e yy according to 
condition (18) and a i: /3i 1. Neglecting these small terms, we obtain 

(n 2 — e xx ) E lxk — e xz E lzk = 0, (40) 
(n 2 - e yy ) E lyk + (-n y n z + e yz ) E lzk = 0, 
£ Z xE\ xk + (—n y n z + e zy ) E lyk + (n 2 — e zz ) E\ zk = 0. 

Below, we will consider the dispersion relation for this system in the collisionless (d 1) 
and collisional (d 3> 1) cases. 
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9.1. Collisionless case 

Consider the case 

(41) 

where an estimation of d is d ~ u ie uj/k 2 c 2 for 5 ~ ctj/mj. Then the dispersion relation has 
two solutions when — a TO = and n 2 — a yy ^ (see the second equation in the system 
(40)). 



The case n 2 z — a yy = This case describes the Alfven waves u 2 = k z c 2 A , where 
ca = B / (Annirrii) 1 ^ 2 . Polarization of the electric field has the form Ei yk ^ and 
Ei x k = E lzk = 0. 



The CclSG 72^, ^yy 

In this case, the dispersion relation reduces to e zz ~ 0. 
Contribution of all other terms can be shown, using expressions (31) and (37), to be small. 
The last equation means that 

5-0. (42) 
Polarization of these perturbations is the following: 

E lzk ^ 0, (43) 

F - Uy F 

n z 

771 _ 771 

(n 2 - a xx ) 

We see that the electric field in the plane of the wave vector is a potential one. However, the 
component E\ xk can be large, E\ xk ~ ctj (k y k z oj C i/k 2 uj) E\ zk . Thus, in general, perturbations 
have an electromagnetic nature. 
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9.2. Collisional case 

In the collisional case, 

d > 1, (44) 
the dispersion relation of the system (40) takes the form 

(n 2 z - a yy ) a zz = 0. (45) 

We see from Equation (45) that the same cases, n 2 z — a yy = and a zz = 0, are possible 
analogously to the case d < 1 considered above. However, polarization in the case 
n z - a yy = ® is different: E lzk = 0, E lyk ^ 0, and E lxk = (n y n z d/a zx ) E lyk ~ {h> ie /u cj ) E lyk . 
This is a mixture of the Alfven and magnetosonic waves. For the dispersion relation a zz = 0, 
the electric field polarization is the following: 

E lzk ± 0, (46) 

IT. _ a xz rp 

(n 2 - a rr ) d 



[n y n z d — a yz ) 
(n 2 z - ciyy) d 



Ei\yk — 7 o VJ~ lzk ' 



where an estimation for d is given above. Thus, the electric field perturbation has all three 
components. 

We see that in the case e zz = the electric field perturbation along the magnetic field 
plays the crucial role in the fast sound speed regime (20). 



10. SOLUTION OF DISPERSION RELATION e zz = 

We now consider the dispersion relation (42) which is suitable also for the collisional 
case. Using notations (23), (32), and (A16)-(A20), we obtain the following dispersion 
relation: 
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( 9 

T e0 I 7 7^ + Q X + ^Te - fir 



d 

+ T i0 ['j— + fi T i - fir. 



d_ 
di 



(47) 



— + f2 x + Vt Te 



2a 



— Qi e (Tio — T e o) 

d 



d_ 

di 



d 



Of 







fine - (7 - 1) 



d_ 

di 



0, 



where d/dt = —iu and fi x = (7 — 1) (Xeo/^eo) k 2 z . All terms fi are defined by the system 
(14). We see that the first four terms on the left hand-side of Equation (47) are symmetric 
concerning the contribution of electrons and ions. The third and fourth terms have appeared 
due to dependence of collision frequency v £ ei ie (nj je o,T e ) on density perturbation. The last 
three terms are connected with perturbation of this frequency because of the electron 
temperature perturbation. All terms proportional to fi e i,i e and fir e j,ie are connected with 
the energy exchange in the thermal equations (3) and (4). We see that taking into account 
the collision frequency perturbation in a general case T e0 7^ T i0 results in considerable 
modification of the dispersion relation. However, this effect should be involved because 
the absence of the thermal equilibrium between electrons and ions can occur, for example, 
in galaxy clusters (e.g., Markevitch et al. 1996; Fox & Loeb 1997; Ettori & Fabian 1998; 
Takizawa 1998). We further use the fact that n i0 = n e0 in our case and v s eiie ~ T e ~ 3 ^ 2 
(Braginskii 1965). Then, we obtain Vl ei = Vl ie and Vl Tei = Vl Tie = — 3fi ie (T e0 — T i0 ) /2T e0 . In 
this case, Equation (47) takes the form 



T e o 7 



d_ 
di 



l^v ~T~ ^T"e — 



d_ 
di 



+ n Ti + 2^ 



(48) 



d 

+ Tio ( 7t^ + fir; - fini 



d_ 
di 



+ fi x + fi Te + 2fii 



+ fiie (^~e0 — Tio) (fi x + fiy e — fi^j 



- fij e (T e Q — Tjo) 



7 | + ^ 



J-eO/ J-eO 



= 0. 
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Below, we will consider different limiting cases of Equation (48). 

10.1. The case Vt ie = 

If we do not take into account the energy exchange, Qi e = 0, and set T e0 = T i0 , then we 
obtain equation 

2 7 w 2 + i [( 7 + 1) (n x + fi Te + n Ti ) - n ne - n ni ] u (49) 

— 2 (Q x + fi^e) &<Ti + ^ne^Ti + fini (fix + ^Te) = 0. 

It is followed from Equation (49) that the ion cooling-heating function modifies the growth 
rate which takes place without this function. Neglecting the contribution of the ion cooling 
and heating, Q Ti = Q ni = 0, we have 

W = -7T" [(7 + 1) («x + ^e) - fine] ■ 
27 

It is easy to see that this solution is a mixture of isochoric and isobaric solutions (Parker 
1953; Field 1965) because we have taken into account the ion temperature perturbation. If 
we neglect the latter, i.e. neglect the second term ~ T i0 in Equation (48), we obtain the 
usual isobaric solution 

U! = (fi x + fire ~ fine) • 

7 

We also see from Equation (48) that for the short-wavelength perturbations when 
Q x ^> u, fire, fine the thermal instability can arise due to the ion cooling function 

oj = — — — ■ — — [(T e0 + T i0 ) Q Ti — T i0 Vt ni ] . 
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10.2. The case Vl ie = oo 

When the frequency f2 ie is much larger than other frequencies, 2Q ie ^> d/dt,Q x ,QTe,i, 
and T e0 = T i0 , then the dispersion relation becomes the following: 



U — — — (Q x + fi^e — ^ne + ^Ti — &"ni) 
27 



(50) 



This is isobaric solution with the electron and ion cooling. 

For the different temperatures of electrons and ions, T e0 7^ T i0 , we obtain 



2 



''io - ( t + ) /; 

J-eO 



u = (3T e0 - T i0 ) (Q x + fl Te ) 



(51) 



In the case T e0 3> T i0 , this equation takes the form, 



UJ = -- [6 (Q x + fiTe) - fine - 5fi Ti ] 



In the opposite case, T e0 <C T i0 , we obtain 



W — (HjPj — fi n j) 

7 



10.3. General case 



In a general case, Equation (48) can be written in the form 



g u + ig x bj -5-2 = 0, 



(52) 



where 



go = 7 (T e0 + T l0 ) , 

gi = [(7T10 + T e0 ) (Q x + ^Te) + (7^0 + ?io) ^Ti - T ef) Vt ne - T i0 fi r . 



1 



T eQ + T iQ 4 + 3 



iO 



eO 



a 



+ (3T e Q — Tjo) a e {Q x + ^Te) — 



^T e0 - 3T J0 (l + ^ 



1 / 
~~ O (TeO + 3Tjo) fiie I fine + ^niT^T - 
^ \ J-eO 



(53) 



Equation (52) can be solved numerically for different cooling functions. 



11. DISCUSSION 

From the results obtained above, we can estimate the relative perturbations of the 
number density and pressure in the fast sound speed regime. Using equations (24)- (27) and 
keeping the main terms, we find expressions for V • Vji, 

—V • Vii ~ — (a,F az + f3 e F elz ) , (54) 
at oz 

V ' Vel ~ ~d~z ( aeFelz + ^ iFil ^ ■ 

Thus, the number density perturbation riji ~ V • v^i is determined by the longitudinal 
electric field E z . If we use the condition of quasineutrality, V • Vji = V • v el , and relation 
rriiFn = — m e F el , we obtain the dispersion relation (42). From equations (21) and (22), 
taking into account polarizations (43) and (46), we have 

dP n dF ilz . . 
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where the value Pji is connected with the pressure perturbation through the relation 

Pn = —L-?£. (56) 
m,jnj at 

From the linear continuity equation and Equations (54)-(56), we find 



d 2 nn 


( OLi 




1 d 2 pn 


dt 2 1 


K m i 


m e j 


1 dz 2 


d 2 n el 


( a e 




I d 2 p el 


dt 2 ~ 1 




rriij 


' dz 2 



It is followed from here the relation between p iX and p ei , 

OLi /3 e \ ( a e (3i 

Pil = J Pel- 



The sum of pressures, p e \ + pn, is equal to 

Pel + Pil = <H Pel < Pel- 

\rrii m e J 

Thus, the total pressure almost does not change. An estimation of the value ctj/mj — f3 e /m e 
gives (see notations (23)) 

oii e \ 1 ( d \ ~ 2 d 2 



mi m e J T \dz J dt 2 
(T e0 ~ T i0 ~ T ). Thus, we obtain nu/n Q ~ pn/po- 

Our dispersion relation (47) does not depend on the wave vector k (except Q x ) and 
magnetic field B . This independence on k is connected with the limiting case (20). The 
absence of the magnetic field is a result of the main contribution to the dispersion relation 
for the condensation mode of the longitudinal perturbed current j\ z under the action of the 
longitudinal electric field E lz (the term e zz Ei z k in the last equation of the system (40)) 
which does not contain the magnetic field. Therefore, a wide spectrum of wavelengths along 
and across the magnetic field can be generated in the framework of conditions (18)- (20) 
and (33). It is worth to note that condition (20) is a particular case which is considered in 
this paper. Other conditions used above are easily satisfied in many astrophysical objects 
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because the ion cyclotron frequency is much larger than the dynamical and ion-electron 
collisional frequencies and the ion Larmor radius is much smaller than the perturbations 
wavelengths. 

The general form of the dispersion relation (48) including the thermal exchange between 
electrons and ions and their different cooling-heating functions and temperatures permits 
us to consider various cases which can be realized in real situations. In particular, Equation 

(49) is available for a weak thermal coupling and equal temperatures, while Equations 

(50) and (51) are appropriate in the case of strong coupling. The intermediate case is 
described by Equation (52), where coefficients (53) contain both different temperatures 
and cooling functions. The large thermal conductivity fl x stabilizes some perturbations. 
We note that one particular Field length (Field 1965) can be obtained from the condition 
Q x = —Vtre + ^ne- However, in the limit fl x ^> uj, Qre,i, ^ne,i, 2il ie , Equation (48) describes 
instability due to the ion cooling (see section 10). 

We have shown that unstable perturbations have an electromagnetic nature (see 
Equations (43) and (46)). Thus, a consideration of potential perturbations given in some 
papers (see section 1) is not adequate. 

From a system of equations (40), we see that for condensation mode e zz = in the 
fast sound speed regime (20) the longitudinal electric field E lz plays the crucial role. The 
transverse wavelengths of unstable perturbations can be both larger and smaller than the 
longitudinal ones. 

The contribution of collisions between electrons and ions in the momentum equations 
depends on the parameter d defined by Equation (38). In both limits (41) (d <C 1) and (44) 
(d ^> 1), the dispersion relation has the same form e zz = 0. 



12. ASTROPHYSICAL IMPLICATIONS 



The fundamental purpose of this paper is to investigate the thermal instability taking 
into account the real multicomponent nature of medium appropriate, e.g., for galaxy 
clusters in a straightforward manner, using the E-approach. Therefore, the different 
cooling-heating functions for the electrons and ions, the energy exchange between them, and 
perturbation of the energy exchange collision frequency for the non-equilibrium background 
state need to be appropriately included. The last effect considered for the first time for 
thermal instabilities contribute to the dispersion relation with the same order of magnitude 
as other effects and considerably modifies it (see, e.g., Equations (50) and (51)). We note 
that the effects mentioned above are not considered in the astrophysical literature using the 
MHD-approach. The key condition for the results obtained is condition (20). It denotes 
that we consider sufncienlty short-wavelengths along the magnetic field. In the simplest 
case without the energy exchange and ion temperature perturbations, the growth rate is in 
the region of the isobaric solution 

LO — (fL, + VLfe — &"ne) 

7 

(see section 10.1). The same growth rate for the condensation mode, we have, for example, 
for the neutral and magnetized media when the perturbation wavelength tends to infinity 
(e.g., Field 1965; Heyvaerts 1974; Loewenstein 1990; Balbus 1991). Thus, from Equation 
(48), we can recover the classical result. 

Conditions on the perturbation scale lengths transverse to the magnetic field are given 
by inequalities (18), (19), and (33) which can be written in the form 
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We see from the first condition that the relation between k± and k z can be arbitrary 
because uj 2 ci ^> u 2 , u ie u, i.e. the cases k± < k z and k± ^> k z are permitted. We emphasize 
that this conclusion is a consequence of the fast sound speed limit (20). The case k z — > 
is not suitable for this limit. The second condition means that the ion sound (or Larmor 
for Q ie = 0) radius is much smaller than the transverse wavelength. For parameters 
appropriate for galaxy clusters, B ~ 1 fiG and T ~ 3 keV, we have u C i ~ 9.6 x 10~ 3 
s _1 and c s ~ (2T/m i ) 1 ^ 2 ~ 7.6 x 10 2 km s _1 (for protons). Thus, the lower limit to the 
transverse wavelength Aj_ = 2n/k± is equal to Ai > 5 x 10 5 km. We note that since the 
growth rates do not depend on k± there is no a specific connection between k±and k z for 
the most unstable modes. The upper limit to k z is determind by the thermal conductivity. 

The dispersion relation (47) for condensation modes does not contain the magnetic field 
because the main dynamics of species in the case under consideration is along the magnetic 
field lines and determined by the longitudinal electric field E z . This field is not captured 
by the ideal magnetic induction equation. We have obtained the dispersion relation which 
cannot followed from the MHD consideration. At the same time, the simplest condensation 
mode growth rate described by Equation (47) coincides with the MHD result (see above 
this section). 

We further shortly outline some important points of our investigation for possible 
observations. We have found the growth rates which contain the separate contribution of the 
cooling functions of electrons and ions in a clear form. It is obvious that both components 
(in multicomponent media, the dust grains and neutrals can also be present) can result in 
the thermal instability in the similar manner. This fact considerably extends possibilities 
for a medium to become unstable. In this connection, it is important to know the functional 
dependence of the cooling functions on the temperature and density for each species. The 
range of the scale lengths of unstable perturbations can enlarge due to contribution of other 
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species except electrons to the instability. For example, the short-wavelength perturbations 
which must be stable because of the large electron thermal conduction can be unstable due 
to contribution of ions to cooling of medium (see section 10). In the fast sound regime when 
the growth rates do not depend on the magnetic field and wavelengths of perturbations, 
both filaments and clouds and pancakes can be observed. It is important to have in mind 
that in this regime the electric field and electric currents of species along the magnetic field 
lines have the crucial role. This result could be used for the diagnostics of the magnetic 
field direction in unstable domains. 



13. CONCLUSION 



We have treated the thermal instability in the electron-ion magnetized plasma which 
is relevant to galaxy clusters, solar corona, and other two-component plasma objects. The 
multicomponent plasma approach have been applied to derive the dispersion relation for 
the condensation modes in the case in which the dynamical frequency is much slower 
than the sound frequency. Our dispersion relation takes into account the electron and 
ion cooling-heating functions, collisions in the momentum equations, energy exchange 
in the thermal equations, different background temperatures of the electrons and ions, 
and perturbation of the energy exchange frequency due to density and temperature 
perturbations. Different limiting cases of the dispersion relation have been considered 
and simple expressions for the growth rates have been obtained. We have shown that 
perturbations have an electromagnetic nature. The important role of the electric field 
perturbation along the background magnetic field has been demonstrated. We have found 
that at conditions under consideration, the condensation must occur along the magnetic field 
lines while the transverse scale sizes can be both larger and smaller than the longitudinal 
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ones. General expressions for the dynamical variables obtained in this paper can be applied 
for a wide range of astrophysical and laboratory plasmas also containing the neutrals and 
dust grains. The results obtained can be useful for interpretation of observations of the 
dense cold regions in astrophysical objects such as IGM, solar corona, and so on. 

I would like gratefully to acknowledge Dr. Mohsen Shadmehri for valuable discussions 
and suggestions and also the anonymous referee, whose useful comments helped in improving 
the manuscript. 



A. APPENDIX 
A.l. Perturbed velocities of species 

In the linear approximation, Equation (1) for the perturbed velocity takes the form 

-7— h b Vji x B , (Al) 

ot mjTijo rrijC 

where pji = n^T^ + rijiTj . From this equation, we can find solutions for the components 
of Vji. For simplicity, we assume that d/dx = because a system is symmetric in the 
transverse direction relative to the z-axis. Then, the x-component of Equation (Al) gives 

dvji x 

= F jlx +u cj v jly , (A2) 

where u C j = qjB /mjC is the cyclotron frequency of species j. Differentiating Equation (Al) 
over t and using Equation (2) in the linear approximation and Equations (15), (16), and 
(A2), we obtain for the y-component of Equation (Al) 

d 2 2 \ dPj! 



-30 - 



where 



Pel 
Pil 



G 2 d 
Dm e dt 

G 3 d 
Dnij dt 



V • Vii + 

V • V e i + 





G 1 d 


m e 


Dm e dt 


Tio 


G 4 d' 


rrii 


Drrii dt / 



V • V e i, 

V-v a . 



(A4) 



The value Pji is connected with the pressure perturbation (see Equation (Al)). Using 
Equations (A2) and (A3), we find 



d 



uj C jdt 



d 2 2 



dPn 
dy 



(A5) 



In Equations (A3) and (A5), notations 



Qjlx ^cjFjly~\ 



dF. 



' ? ' dt ' 

dFj lx 



(A6) 



dt 



are introduced. We see from these equations that the thermal pressure effect on the velocity 
v ilx is much larger than that on v ily when d/dt <C u ci . The ^-component of Equation (Al) 
can be written in the form 

dP» + dFju ^ 



dt 2 



dz 



dt 



A. 2. Calculation of V • Vji and Pji 

We have 

dvjiy dv jlz 



V-v, = ^ + ^f. (AS, 



Using Equations (A3), (A4), (A7) and (A8), we obtain 



Lu V • V e i + L 2e V • V a = if e i, 

LuV ■ Vn + L 2i V ■ v e i = # a . 



(A9) 
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Here, 



d 3 Q, 



dydt 2 



+ 



d 2 

W 2 + ^ 



d 2 F n 



JlZ 



dzdt 



2 \d^_ 

dt 2 +UJ ce) df2 



2 1 " r 



and operators L ± j and L 2 j are the following: 

Lie 

Lu 

L>2e 
L 3 j 



d 2 M 2 \d 2 (T t0 



G 1 d 
Dm e dt 
G A d 



j ° 2 9 L 

J^3e— 77T> ^2i 



rrii Drrii dt 
G 3 d 



Dm e dt ' 01 Drrii dt ' 

dy 2 dt 2 + \dt 2 +UJcj ) dz 2 ' 



From a system of equations (A9), we find 



(A10) 



(All) 



LV ■ v el = — L 2e Hii + L u H e i, 
LV • Vji = —L 2 iH el + Li e 



(A12) 



where 



L — L le L u — L 2e L 



2i- 



(A13) 



The values Pji and P ei can be found, substituting solutions (A12) into expressions (A4), 

G 4 d 



' G:; d ^ /7a, 

Dwij dt e \m» Umj <9t 

G 3 8^ /Tip g 4 9 

Dm; <9t * \mj Z)mj dt 



L 2 i 



Ha 

H e li 



(A14) 



LPel = 



' G 2 d_^ | /T e0 
Dro e dt e \ m e 



d d 
Dm e dt 
d d 
Dm e dt 



Lu 



■>2e 



H e l 
Hi\. 



(A15) 
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A. 3. Expressions for D and Gi,2,3,4 



We now give expressions for values given by a system (17): 
D = 



d \ ( d \ d 2 ( d \ d 2 

1 Q x + ^Te) ( 1 - + n Ti + n ie ) — + (n ei + n Tei ) [— + n Ti 



dt 



dt 



dt 2 



dt 



dt 2 ' 



(A16) 



G l =T ( 



eO 



d 

n ne - (7 - 1) — 



^ + fi Tl + n ie ) ^ + n ei (T e0 - T, ) + ^ |r, (A17) 



dt 



G 2 — Vt ei TiQ 



tt ni - (7 - 1) 



d_ 

dt 

d 



tt ne - (7 - 1) — 



^ + fi« (T e0 - T l0 ) (JL + ^ |, (A18) 
- - a e (T e0 - T l0 ) ( - + fi x + fi Te J -, (A19) 



d_ 

dt 



d 

tt ni - (7 - 1) — 



. 9 \ 9 

— fije (T e0 — Tj ) ( — + £\ + fij-e J 7^- 



(A20) 



A. 4. Simplification of Equations (A14) and (A15) 

We shall further calculate coefficients by Hji in Equations (A14) and (A15). Using 
expressions (All), we find 

G 3 d fT i0 G 4 d\ _ G 3 (d 2 2 \ d 3 
D^d~t U + {m t ~^m) L2l -D^ t {w +Uci JW {A21) 

and 

«3 e + _ ft a n 1 /iv, _ cu a w* + x * + 1 



Dm, : <9t \ Dm, : <9t / D \ rrii m, : dt J \ dt 2 ce J dt 2 Dm e rrii 



(A22) 

In Equation (A22), we have introduced notation 

K = ^ (G 2 G 3 - G 1 G 4 ) + (T e0 G 4 + T i0 Gi) ^ - DT e0 T l0 . (A23) 
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Calculations show that the value (G 2 G 3 — G1G4) has a simple form, i.e. 



— (G2G3 — G1G4) — Vlie (T e0 — Tj ) T f 



ef) 



d 

tt ne - (7 - 1) — 



+ fi e j (TiQ — T e0 ) T i0 



d 

tt ni - (7 - 1) — 
(A24) 



— TeoTiQ 



fine - (7 - i; 



d_ 
dt 



d 

tt ni - (7 - 1) — 



Using expressions (A16), (A17), (A20), and (A24), we obtain for the operator K (A23) the 
simple form, 



K = -n ie TlW e ^ - (fi ei T l0 + n Tei T e0 ) T i0 W~ - T e0 T i0 W e W~, (A25) 



where notations 



d_ 

dt 
d_ 

dt 



a 

W e = 7T^ + fi X + fi^e - fine, 



Wi = 7^7 + fiTi - fi, 



(A26) 



are 



introduced. Using Equations (A21) and (A22), Equation (A14) for Pn takes the form 



DLP a 
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Analogous consideration of Equation (A15) leads to the following equation for P el : 



DLP, 



el 
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Operators 
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jjTjQ ^ G4 d 

mi rrii dt 



(A27) 



(A28) 
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can be found by using Equations (A16), (A17), (A20), and (A26) 

D —--i = ~ W * (I + ^ + n*) ti (A29) 
m e m e at m e \ot J at 2 

1 (d \ d 2 



D T M G,d_T i0 w fd ,0,0 V 2 



A. 5. Operator L in a general form 



Using expressions (All), we find from Equation (A13) 



where 



M-N 1 L 3e L 3i K, (A30) 



m e niiD 



M=i^+<^)(^+<^)^z, (A31) 



dt 2 ce J \dt 2 cl J dt v 
v , 9 2 2 \ d 2 fT e0 G 1 d\ (d 2 2 \d 2 (T l0 G 4 d' 
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